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Atoms trapped in magic wavelength optical lattices provide a Doppler- and collision-free dense 
ensemble of quantum emitters ideal for fast high precision spectroscopy and thus they are the basis 
of the best optical clock setups to date. Despite the minute optical dipole moments the inherent 
long range dipole-dipole interactions in such lattices at high densities generate measurable line 
shifts, increased dephasing and modified decay rates. We show that these effects can be resonantly 
enhanced or suppressed depending on lattice constant, geometry and excitation procedure. While 
these interactions generally limit the accuracy and precision of Ramsey spectroscopy, under optimal 
conditions collective effects can be exploited to yield zero effective shifts and long dipole lifetimes 
for a measurement precision beyond a noninteracting ensemble. In particular, 2D lattices with a 
lattice constant below the optical wavelength feature an almost ideal performance. 


Since the turn of the century the technology of manip¬ 
ulating and controlling ultracold atoms and molecules 
with laser light has seen breathtaking advances [1-3]. 
Following the seminal first demonstration of a quantum 
phase transition in an optical lattice [4], nowadays the 
so-called Mott insulator phase for atoms in an optical 
lattice can be prepared almost routinely [5, 6] and exper¬ 
iments with photo-associated ultracold molecules have 
reached a comparable level of control [7 -10]. Employ¬ 
ing the particles’ internal structure coherent interactions 
between the atoms at neighboring sites in such a lattice 
can be tailored to a large extend, e.g. via spin-dependent 
tunneling [11]. 

In one of the typical experimental setups atoms pos¬ 
sessing a long-lived clock transition are prepared in an 
optical lattice using a differential light shift free (magic) 
trapping wavelength [12, 13]. Its most prominent appli¬ 
cation is the implementation of the world’s best optical 
clocks [14-16]. 

When the atoms in such a lattice are excited on an 
optical or infrared transition they will interact stronger 
and on a much longer range via dipole-dipole energy ex¬ 
change than via tunneling or collisions. At sufficient den¬ 
sities the dipole interaction strength surpasses the ex¬ 
cited state lifetime and resonantly enhanced collective 
excitations analogous to excitons in solid state physics 
appear [17, 18]. For polar molecules in optical lattices, 
such long wavelength collective excitations even domi¬ 
nate the dynamics [19] and can form the basis for study¬ 
ing generic phenomena of solid state physics by means 
of tailored toy models [1]. In the case of clock transi¬ 
tions, the extremely tiny dipole moment keeps these in¬ 
teractions small in absolute magnitude. Nevertheless, for 
densities close to unit filling the exciton’s effective transi¬ 
tion frequencies and their spontaneous decay is governed 
by dipole-dipole interaction [20] and deviates from the 
bare atom case. 

In a standard Ramsey interrogation sequence used as 
the generic basis for a clock setup, the first pulse creates 
a product state of all spins prepared in a 50% super¬ 
position of ground and excited state with equal phase 



FIG. 1. (Color online) Scheme of a 2D optical lattice filled 
with clock atoms interacting via dipole-dipole energy ex¬ 
change flij and a collectively modified spontaneous emission 
r ij. In a mean field treatment these quantities can be summed 
up yielding the two effective quantities fi eff and F Gff , which 
suffice to approximate the system dynamics. 


and with all dipoles in parallel alignment. This state 
features the maximally possible dipole moment and thus 
exhibits superradiance. Even a tiny single particle spon¬ 
taneous emission rate can be strongly enhanced, so that 
collective decay will become the dominant factor limit¬ 
ing measurement time and precision for a large and dense 
ensemble [21]. While in the current setups based on ID 
lattices with fairly low filling this perturbation is largely 
negligible [22], we will see that close to unit filling, it be¬ 
comes sizable and constitutes a dominant effect in 2D and 
3D lattices. As a key quantity for characterizing the col¬ 
lective dynamic, we calculate the decay and phase shift 
of the collective dipole between the two Ramsey pulses, 
which determines the contrast of the Ramsey sequence. 
Note, that due to the pairwise nature of the interaction, 
rephasing pulses cannot correct these errors easily. 

For a consistent quantum treatment of the open sys¬ 
tem dynamics, one has to solve the corresponding mas¬ 
ter equation for the atomic density matrix p obtained by 
tracing over the electromagnetic vacuum modes [17, 23], 
which reads 


p = i[p,H] + C[p). 


(1) 
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As the particle Hilbert space grows exponentially with 
their number, only small ensembles far from any realistic 
size can be treated rigorously. However, since for preci¬ 
sion measurements we need to evaluate collective effects 
to many digits, reliable and converging numerical mod¬ 
els are required to obtain meaningful predictions. For 
rather small atom numbers (< 12) a numerical solution 
of the full master equation yields non-negligible shifts al¬ 
ready [20]. For larger ensembles at very low densities, a 
cluster approach based on statistical averaging over small 
particle number configurations has produced first esti¬ 
mates of the scaling of the dephasing with the system’s 
size and density [14]. Recently, important self synchro¬ 
nization effects through dipole coupling were studied in 
a very high density limit using simplifying assumptions 
for the coupling [24]. 

In this letter we present an extensive analysis of the 
collective dynamics for fully populated regular lattices of 
different geometries in one to three dimensions containing 
a large or even an infinite number of particles. Our cen¬ 
tral goal is to estimate the magnitude of the dipole phase 
shift and collective decay as a function of the lattice and 
excitation geometry. Surprisingly, besides a resonant en¬ 
hancement of decay and clephasing, we find cases where 
collective effects lead to improvements of the expected 
measurement precision beyond the independent particle 
level by virtue of subradiant excitonic states. Numeri¬ 
cally we start from a mean field approach whose validity 
has been extensively tested in recent work [25]. Due to its 
simplicity we can apply it to very large particle numbers 
well approximating even infinite size systems. Its accu¬ 
racy, however, breaks down at very close distances, as 
it, obviously, cannot capture the full correlations among 
all particles. Similar deliberations in the classical regime 
have recently been put forward [26]. 

Model ~ We consider an ensemble of N identical ef¬ 
fective two-level atoms with transition frequency wo and 
lifetime 7 at positions 77 ( i = 1..N) interacting via dipole- 
dipole coupling described by the Hamiltonian 

H = *7- ( 2 ) 

Here, af denotes the raising (lowering) operator of the 
i-th atom and = 17 G (fcoA) represents the coherent 
energy exchange among the atoms with ko = ujq/c = 
27t/Ao and ry,- = |rj — rj\ being the distance between 
atom i and j. Collective spontaneous emission is ac¬ 
counted for by a Lindblad-type Liouvillian of the form 

c Ip\ = 2 zL r y (Ci)( 2<j~pa+ - af a~p - pa+ a ~), (3) 

i,j 

where the off-diagonal spontaneous emission rates F^ = 
§ 7 F (koVij) introduce super- and subradiant decay chan¬ 


nels [27]. Explicitly we have 

sin ^ /cost; sin A 

F{0 = a-^+L3[-jr--gr) (4a) 

cos£ n /sin£ cos A . 

G(0 = —a—Gf + /3 + yj (4b) 

with a = 1 — cos 2 9 and ,0 = 1 — 3 cos 2 0, where 9 repre¬ 
sents the angle between the line connecting atom i and 
j and the atomic dipole orientation, which is assumed 
equal for each atom. 

Mean field approximation - To enable a study of large 
particle numbers we derive the equations of motion for 
the expectation values of the Pauli operators for the k- 
th atom as detailed in the supplement [28]. Assuming a 
separable density operator, these equations can be trans¬ 
formed to a closed form by factorizing the two-particle 
correlations via (cdVp « (crf)(crj) for p, v £ {x,y,z}. 
As shown in a previous publication [25] these reduced 
equations still capture the major part of the interaction 
up to a moderate interaction strength. For more accu¬ 
rate studies we can add second order corrections, which 
increase the computational effort but are still tractable 
and allow for estimating the validity of the mean field 
approximation. 

Symmetric configurations - In the case of special sym¬ 
metric geometries with each atom initially in the same 
state and subject to the same effective interactions with 
the other atoms, the equations of motion become identi¬ 
cal. It then suffices to solve the following relatively simple 
system, 

(a x ) = ft ef V)<A> - \ (7 - r ef V>) <^>, (5a) 

<A> = -n ef V)<A) - \ (7 - (A>) (a v ), (5b) 

(A) = -7(1 + (A)) - V' ((A ) 2 + (A) 2 ) . (5c) 

With this, the numerical difficulty reduces from solving 
a huge set of coupled nonlinear equations to calculating 
the effective couplings for a large number of identically 
evolving particles, i.e. 

N N 

f r s = ^n r eff = E r H’ (6) 

3 =2 3 =2 

Of course, a rigorous symmetry condition is fulfilled for 
very few atomic distributions only, e.g. in regular poly¬ 
gons and cuboids but also and probably with more rel¬ 
evance for actual experimental setups, in infinite chains 
and infinite periodic 2D and 3D lattice configurations. 
The essence of the interactions within the entire lattice 
structure is captured solely by two real numbers, the ef¬ 
fective coupling fl cff and the collective decay rate F eff . In 
an atomic clock setup one seeks to minimize the energy 
shifts fl cft as far as possible and find configurations with 
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a maximally negative r eff , thereby minimizing decay and 
allowing for an as long as possible interrogation time. 

Finite systems - Let us now turn to concrete examples 
for which the effective quantities can be calculated easily. 
The most obvious structures where all particles exhibit 
the same dynamics is any regular polygon. This might 
not be the most practical example but displays the un¬ 
derlying physical mechanisms nicely. In Fig. 2 we com¬ 
pare the effective interactions for a square, a ten-sided 
and a 100000-sided polygon. The square (red) shows a 
behavior quite similar to the underlying functions F(tf) 
and G(£), while the two larger polygons (blue and green) 
exemplify a more complex scenario. Note their strong 
size dependent variation and the divergence of the ef¬ 
fective interaction strengths at integral values of dAo for 
large numbers that emerges due to the accumulation of 
l/£ contributions. We see that even with relatively large 
atom spacing on the scale of the transition wavelength, 
cooperative collective effects are large and vary strongly 
with distance. 



0 1 2 30 1 2 3 

d/ Ao d/ Ao 

FIG. 2. (Color online) Distance dependence of the effective 
dipole coupling fi ofi and r cft for a square (red), a ten-sided 
(blue) and a 100000-sided (green) regular polygon. The fewer 
particles the closer the functions resemble the underlying cou¬ 
plings fi ij and r,;,. The divergences at integral d/ Ao result 
from the l/£ terms in F(£) and G(£). 

Infinite systems - Of course, extended regular sys¬ 
tems, i.e. large periodic optical lattices, are experimen¬ 
tally more relevant. At this point our mean field treat¬ 
ment opens up the possibility to investigate realistic atom 
numbers numerically. Fig. 3 depicts these effective cou¬ 
plings as a function of the lattice constant for an infi¬ 
nite chain, a square lattice and a hexagonal lattice. For 
comparison, we have added the results for smaller atom 
numbers to demonstrate cutoff effects, which might al¬ 
low for an estimate of the necessary particle numbers 
for actual experimental setups as well. Again, we ob¬ 
serve strong variations and divergences of the couplings 
at integral values of d/Ao that stem from the long range 
l/£-terms in the functions F(£) and G(£). These mani¬ 
fest themselves in a much more pronounced way at huge 
atom numbers and therefore underpin the importance of 
properly including long range couplings. 

Note that for the two-dimensional square lattice and 


the hexagonal lattice r cff exhibits a broad minimum for 
the effective decay close to F eff = —1 for lattice constants 
d < Ao- In this case atomic decay is strongly inhibited 
by the neighboring atoms. This immediately suggests 
that such two-dimensional setups are very favorable for 
building future generation lattice clocks. For such large 
negative r eff subradiant decay will dominate the system 
dynamics allowing for much longer Ramsey delay times 
and thus offering a higher overall precision [21]. Sim¬ 
ilarly we can identify lattice constants with zero effec¬ 
tive shifts increasing the possible clock accuracy. Upon 
extending these calculations to three dimensional lat¬ 
tices, we find that the necessary atom numbers to obtain 
smooth converging behavior analogous to the ID and 2D 
case are beyond our current numerical capabilities. For 
tractable particle numbers of about 10 12 the resulting 
effective quantities still fluctuate strongly, predicting po¬ 
tential problems for such 3D clock setups. 

Tailoring atomic excitations - So far we have assumed 
a phase-symmetric excitation of all atoms in the first 
Ramsey pulse. In a practical plane wave excitation 
scheme this corresponds to an illumination from the side 
at a right angle. In general, however, the effective cou¬ 
plings D eff and T cff will change, when we allow for a lo¬ 
cal phase shift imprinted on the atoms by the excitation 
laser, which comes into play when they are addressed 
at a different angle. In a tt/2 Ramsey spectroscopy se¬ 
quence [29], where the atoms are initialized in a coherent 
superposition of ground and excited state, the excitation 
phase manifests itself on the excited state directly, i.e. 

N . 

\*) = <£)- / =(\g)+e iA ‘M- 1) \e)). (7) 

j =i v 

In our mean field treatment we can make use of the sys¬ 
tem’s translation symmetry and absorb this phase into 
the effective couplings as discussed in detail in the sup¬ 
plement [28]. For A(f> = 0 we reproduce the above case. 
The closer the phase gets to A <p = ir, however, the more 
half-integral values of dXo seem to yield minimal shifts 
and the maximally negative r off as seen in Fig. 4. This 
can be attributed to the fact that, had the emitted light 
interfered constructively at integral and destructively at 
half-integral distances for A0 = 0, it will do exactly the 
opposite at A<ft = 7r. Furthermore, this result indicates 
that addressing atoms transversally (Af = 0) seems to 
be more favorable at typical magic wavelength trapping 
distances, e.g. d/X q ~ 0.58 for 87 Sr [12, 20, 30]. Let us 
emphasize again, that for d <C Ao the mean field approach 
breaks down. At such small distances one would rather 
turn to the Dicke model [31], which allows for reducing N 
two-level emitters to one effective spin 7V/2-system[24]. 

Now, having obtained effective interaction parameters 
for a wide range of configurations and system sizes, let us 
finally discuss their consequences in typical cases. Fig. 5 
shows the time evolution of an infinite chain that is ini- 




















4 



FIG. 3. (Color online) Distance dependence of the effective quantities fi cff and r efi for an infinite equidistant chain, a square 
lattice and a hexagonal lattice (dashed black) compared to their not yet converged finite counterparts of 10, 4 * 10 4 and 10° 
particles respectively (solid red). Again, we find divergences at integral d /Ao owing to the l/£-terms in F{Q and G(£)- In the 
2D configurations F eff plateaus at —1 for d < Ao, suggesting that this parameter range will be the most favorable for clock 
setups as decay is strongly suppressed. 
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FIG. 4. (Color online) Effective interactions fl eff and r e£I for 
an infinite chain with spacing a where the system is initially 
prepared with all spins in a Ramsey state but with phase shift 
of A <j> between neighboring spins. The dashed lines indicate 
where fi off = 0, which is optimal for an optical clock. 

tialized in a symmetric Ramsey state by a quick n/2- 
pulse, i.e. every atom is in a coherent superposition of 
ground and excited state, with the excitations baring ei¬ 


ther no phase corresponding to a transverse addressing 
of the atoms, or a phase shift of Ac/) = ir between neigh¬ 
boring atoms. The lattice constants have been chosen to 
be approximately Xq/2 as would be typical in an optical 
lattice [5]. Note, however, that choosing exactly A 0 /2 
would give us the l/£ divergence discussed above at ev¬ 
ery second lattice sight, which would, of course, be detri¬ 
mental to useful numerical results. We observe that the 
dipole’s lifetime varies strongly, comparing the subradi¬ 
ant behavior (red) where the collective dipole lives much 
longer than the natural lifetime of the atom to the su¬ 
perradiant (green) regime where the excitation vanishes 
very quickly. Additionally, to highlight the validity of 
the mean held approach, we provide the results of a sec¬ 
ond order cumulant expansion simulation, which demon¬ 
strate a satisfactory agreement confirming the mean held 
approximation. 

Conclusions - We have shown that in densely filled 
magic wavelength optical lattices dipole-dipole interac¬ 
tion and collective decay signihcantly changes the evo¬ 
lution of an induced collective dipole. Due to the long 
range-nature of the coupling, shifts appear even for very 
long lived clock states despite their minute dipole mo¬ 
ment. This has direct consequences on the accuracy and 
precision of a Ramsey interrogation procedure as used 
in state of the art atomic clocks. Shifts and dephasing 
in extended systems strongly depend on dimensionality 
and geometry of the lattice, exhibiting resonant enhance¬ 
ments at particular lattice constants. While at current 
operating densities and numbers these shifts can be ne- 
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FIG. 5. (Color online) Three different examples for the time 
evolution of the Pauli expectation values for a chain with spac¬ 
ing d where initially all spins are prepared in a coherent super¬ 
position of ground and excited state, but have a phase shift of 
A <f>. The parameters used are d = 0.792Ao with A(p = 0 (red 
triangles), where fl cff = 0 and T off is nearly optimal, as well 
as d = 0.49Ao (green squares) and d = 0.51Ao (blue circles), 
both with A 4> = 7r which are close to a P cff discontinuity. The 
solid lines correspond to a solution of a second order cumulant 
expansion model with 200 particles and demonstrate a very 
good qualitative and sufficient quantitative agreement with 
the infinite mean field description. 


glected compared to other technical imperfections, they 
constitute inherent, fundamental perturbations. 

As our central result we have identified optimal opera¬ 
tion geometries, which combine negligible effective shifts 
and a strong suppression of decay, superior to indepen¬ 
dent atom ensembles. In particular, for a ID lattice with 
a tailored excitation angle and for a 2D hexagonal lattice 
optimal operation parameters for a next generation clock 
setup have been demonstrated. In this sense it seems 
possible to implement a high density dark exciton based 
atomic clock geometry with shifts many orders of magni¬ 
tude below a single Hz and almost unlimited exciton life 
times. 

While for most considerations we have focused on the 
case of clock transitions, the same physics is present 
in a more prominent and experimentally easier observ¬ 
able form for broader transitions. Optimizing geome¬ 
tries will also be relevant for devices such as superradiant 
lasers [32, 33] or lattice based optical memories. 
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Atoms deeply trapped in magic wavelength optical lattices provide a Doppler- and collision-free 
dense ensemble of quantum emitters ideal for high precision spectroscopy. Thus, they are the basis 
of some of the best optical clock setups to date. However, despite their minute optical dipole 
moments the inherent long range dipole-dipole interactions in such lattices generate line shifts, 
dephasing and modified decay. We show that in a perfectly filled lattice these effects are resonantly 
enhanced depending on lattice constant, lattice geometry and excitation scheme inducing clock 
shifts of many atomic linewidths and reducing measurement precision via superradiance. However, 
under optimal conditions collective effects can be exploited to yield zero effective shifts and prolong 
dipole lifetimes beyond the single atom decay. In particular we identify 2D hexagonal or square 
lattices with lattice constants below the optical wavelength as most promising configurations for an 
accuracy and precision well below the independent ensemble limit. This geometry should also be 
an ideal basis for related applications such as superradiant lasers, precision magnetometry or long 
lived quantum memories. 


Since the turn of the century the technology of manip¬ 
ulating and controlling ultracold atoms with lasers has 
seen breathtaking advances [1-3]. Following the seminal 
first demonstration of a quantum phase transition in an 
optical lattice [4], nowadays the so-called Mott insulator 
state can be prepared routinely [5, 6]. Experiments with 
photo-associated ultracold molecules have reached a com¬ 
parable control [7-10] and coherent interactions between 
the atoms at neighboring sites can be tailored [11]. 

For some of the world’s best optical clocks [12-14], 
atoms with a long-lived clock transition are prepared 
in an optical lattice using a differential light shift free 
(magic) trapping wavelength [15, 16]. In principle, this 
provides for a Doppler and collision free dense ensem¬ 
ble with negligible inhomogeneous broadening. However, 
when excited optically emitters will nevertheless interact 
via long range resonant dipole-dipole coupling [17]. 

Here we show that at sufficient densities the dipole in¬ 
teraction strength surpasses the excited state lifetime and 
collective excitations analogous to excitons appear [18]. 
For polar molecules in optical lattices they even dominate 
the dynamics [19] and allow for studying generic phenom¬ 
ena of solid state physics [1]. For clock transitions the 
extremely tiny dipole moment keeps these interactions 
small in absolute magnitude. Still, the exciton’s effec¬ 
tive transition frequencies and their spontaneous decay is 
governed by dipole-dipole interaction [20] deviating from 
the bare atom case. This limits accuracy and precision 
of corresponding clock setups. In an idealized Ramsey 
sequence for a clock setup, the first laser pulse creates a 
product state of all atoms prepared in a 50% superpo¬ 
sition of ground and excited state with equal phase and 
all dipoles aligned in parallel. This state features the 
maximally possible dipole moment and typically exhibits 
superradiance. Even a tiny single particle spontaneous 
emission rate can be that strongly enhanced, that collec¬ 
tive decay becomes a dominant factor limiting measure¬ 
ment time and precision [21]. In current setups based 
on ID lattices with low filling this perturbation is often 
negligible compared to other noise like collisions, black 



FIG. 1. (Color online) Scheme of a 2D optical lattice filled 
with clock atoms interacting via dipole-dipole energy ex¬ 
change flij and a collectively modified spontaneous emission 
Tij at two different lattice constants shown in blue and yel¬ 
low. In a mean field treatment with translation invariance the 
sum over all interaction terms yields two effective couplings 
D cff and F eff only, which govern the approximate system dy¬ 
namics. 


body shifts or reference cavity fluctuations. However, in 
lattices with unit filling, dipole-dipole interaction shifts 
are larger than the atomic linewidth and constitute a sig¬ 
nificant inherent perturbation. Note that their absolute 
magnitude scales with the atomic dipole moment and 
thus strongly depends on the chosen transition. 

As a key quantity to capture the collective modifi¬ 
cations of the system dynamics, we use the decay and 
phase shift of the collective dipole generated by the first 
Ramsey pulse, which determines the contrast and shift 
of the central Ramsey fringes. Note, that due to the 
pairwise nature of dipole-dipole interactions a rephasing 
pulse cannot correct these errors. Here we ignore inter¬ 
action induced perturbations during the Ramsey pulses, 
which introduce extra noise but could be reduced by very 
fast pulses or improved sequences [22]. 

For this we numerically solve the well established mas¬ 
ter equation for the atomic density matrix p including 
optical dipole-dipole interaction obtained by tracing over 
















2 


the electromagnetic vacuum field [17, 23, 24], 

p = i\p,H\+C[p\. (1) 

As previously shown for small atom numbers (N < 12 ) a 
numerical solution of the full master equation yields non- 
negligible shifts already [20, 25]. As the Hilbert space 
grows exponentially with atom number, the full equation 
cannot be solved for ensembles of a realistic size. Since 
for precision measurements we need to evaluate collec¬ 
tive effects precisely, reliable and converging alternative 
numerical methods are required. For larger ensembles at 
low densities a cluster approach has produced first esti¬ 
mates of the scaling of the dephasing with the system’s 
size and density [12]. Recently, important self synchro¬ 
nization effects through dipole coupling were studied in 
a very high density limit using simplifying assumptions 
for the coupling [26]. 

In this letter we present an extensive analysis of the 
collective dynamics for fully populated lattices of differ¬ 
ent geometries and sizes containing a large or even an 
infinite number of particles. Our primary goal is to es¬ 
timate the magnitude of the dipole phase shift and col¬ 
lective decay as a function of lattice and excitation ge¬ 
ometry. Besides a resonant enhancement of shifts, decay 
and dephasing, we find cases where collective effects lead 
to improvements of the maximally achievable measure¬ 
ment precision beyond the independent particle level by 
virtue of subradiant states. We concentrate on an ideal¬ 
ized setup ignoring lattice shifts, thermal effects or the 
hopping of atoms. 

Numerically we apply an enhanced mean field ap¬ 
proach whose validity has been extensively tested in re¬ 
cent work [27], which we can scale to realistic numbers 
up to N ss 10 5 particles. If the particle distribution ex¬ 
hibits symmetries numbers up to even 10 10 are possible, 
well approximating infinite systems in ID and 2D. Its 
accuracy, however, breaks down at very close distances 
as it cannot capture higher order correlations. Similar 
deliberations for classical dipoles have recently been put 
forward [28]. 

Model - We consider an ensemble of N identical ef¬ 
fective two-level atoms with transition frequency wo and 
inverse lifetime 7 at positions 77 (i = 1 ..IV) interacting via 
optical dipole-dipole coupling described by the Hamilto¬ 
nian [17, 24] 

H = ( 2 ) 

Here, af denotes the raising (lowering) operator of the 
z-th atom and fly = (korij) represents the energy 
exchange with /c 0 = ojq/c = and ry,- = 177 — 77 ■] 

being the distance between atoms i and j. Collective 
spontaneous emission is accounted for by a Liouvillian of 
the form [17, 23] 

c [p\ = \ T ij( r ij)( 2a i P a t ~ a t a j P ~ P a t a j )> ( 3 ) 

i,3 


where the off-diagonal rates Fy = § 7 F (fcoDj) introduce 
super- and subradiant decay [24]. Explicitly we have 

sin ^ . /cos£ sin A 

= (4a) 

cos£ _ /sin£ cos A , 

G(i) = -a-^ + 0 + yj (4b) 

with a = 1 — cos 2 6 and /? = 1 — 3 cos 2 9 , where 9 repre¬ 
sents the angle between the line connecting atoms i and 
j and the common atomic dipole orientation. 

Mean field approximation - To study large particle 
numbers we derive the equations of motion for the ex¬ 
pectation values of the Pauli operators for the k-th atom 
as detailed in the appendix A. Assuming a separable den¬ 
sity operator and factorizing the two-particle correlations 
via (erferf) « (erf)(crf) for /z, v £ { 2 :, y, 2 } they transform 
to a closed set. As shown previously [27] these equations 
still capture the major part of the interaction up to a 
moderate interaction strength. For more accurate stud¬ 
ies one may add second order corrections increasing the 
computational effort. 

Symmetric configurations - For symmetric geometries 
with each atom initially in the same state and subject to 
the same effective interactions, the equations of motion 
for all particles become identical and read 

(A) = S } eff (A)(A> - \ (7 - r ef V>) <0, (5a) 

(A) = -n^vx**) - \ (7 - r eff (A)) (A), (5b) 
(A) = -7(1 + (A)) - V'((A> 2 + (A) 2 ). ( 5 c) 

Hence instead of solving a huge set of coupled nonlinear 
equations we need to determine the effective couplings, 
i.e. 

N N 

fi eff = 5> y r eff = ^r lj , (6) 

3 =2 j—2 

Of course, such a rigorous symmetry condition is fulfilled 
for very few atomic distributions only. Then, the essence 
of the interactions within the entire lattice is captured 
solely by two real numbers, the effective coupling D ofi and 
the collective decay rate F eff . In a clock setup one seeks 
to minimize the energy shifts and find configurations 
with a maximally negative r eff , minimizing decay and 
allowing for an as long as possible interrogation time. 

Finite systems - Firstly, for finite symmetric configura¬ 
tions the effective quantities can be calculated easily. The 
most obvious symmetric structures are regular polygons. 
This might not be the most practical example but nicely 
displays the underlying physics. In Fig. 2 we compare the 
parameters for a square, a ten-sided and a 100000 -sided 
polygon. The square shows a behavior quite similar to 
the underlying functions F(£) and G(£), while the two 
larger polygons exhibit strong size dependent variations, 
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particularly at integral values of dAo emerging from the 
accumulation of many l/£ contributions. Note that even 
with a relatively large atom spacing, cooperative collec¬ 
tive effects are sizable and vary strongly with distance. 




FIG. 2. (Color online) Distance dependence of the effective 
dipole coupling fl cfI and r cfi for a square (red), a ten-sided 
(blue) and a 100000-sided (green) regular polygon. The fewer 
particles the closer the functions resemble the underlying cou¬ 
plings f lij and Tij. The divergences at integral d/ Ao result 
from the l/£ terms in F(£) and G(£). 

Infinite systems - In practice, extended regular sys¬ 
tems, i.e. large periodic lattices, are experimentally more 
relevant. Fig. 3 depicts the effective couplings for an 
infinite chain, a square lattice and a hexagonal lattice. 
For comparison, we have overlaid the results for smaller 
atom numbers to demonstrate finite size effects. We ob¬ 
serve stronger variations and again divergences at inte¬ 
gral values of d/Ao . These manifest themselves in a much 
more pronounced way at huge atom numbers and there¬ 
fore underpin the importance of properly treating long 
range interactions. 

Note that for the two-dimensional square lattice and 
the hexagonal lattice r eff exhibits a broad minimum for 
the effective decay close to r off = — 1 for d < Ao, where 
atomic decay is strongly inhibited. This favors such two- 
dimensional setups for lattice clocks as subradiant decay 
will dominate the system dynamics allowing for much 
longer Ramsey delay times and thus offering a higher 
overall precision [21]. Similarly we can identify lattice 
constants with a zero effective shift increasing clock ac¬ 
curacy. Extending these calculations to three dimen¬ 
sional lattices, we find that the necessary atom numbers 
to obtain smooth converging behavior are beyond our 
current numerical capabilities. For particle numbers of 
about 10 12 the resulting effective quantities still fluctu¬ 
ate strongly, predicting potential problems for such 3D 
clock setups. A demonstration of this effect can be found 
in the appendix C. 

Tailoring atomic excitations - So far we assumed a 
phase-symmetric excitation of all atoms by the first Ram¬ 
sey pulse. In a practical excitation scheme this corre¬ 
sponds to illumination at right angle. In general, how¬ 
ever, the effective couplings D cff and r eff will change, 
when we allow for a local phase shift imprinted on the 
atoms. In a tt/2 Ramsey sequence [29] the excitation 


phase appears on the excited state directly, i.e. 
n 1 

\*)=®-^(\9)+e iA W-V\e)). ( 7 ) 

j —i v 

In our treatment we can exploit the system’s symmetry 
and absorb this phase into the effective couplings (Ap¬ 
pendix B). For A (f> = 0 we recover the above results. 
The closer the phase shift gets to A (f> = n, however, the 
more half-integral values of d\o yield minimal shifts and 
the maximally negative F e ® as seen in Fig. 4. Since the 
emitted light has interfered constructively at integral and 
destructively at half-integral distances for A (f> = 0, it will 
do exactly the opposite at A = it. Furthermore, ad¬ 
dressing atoms transversally (A (f> = 0) seems more favor¬ 
able at typical magic wavelength trapping distances, e.g. 
d/ Ao ~ 0.58 for 87 Sr [15, 20, 30]. Again, for d <C Ao the 
mean field approach breaks down and one should rather 
turn to the Dicke model [31], reducing N two-level emit¬ 
ters to one effective spin A/2-system [26]. Let us finally 
discuss the consequences for typical cases. Fig. 5 shows 
the time evolution of the average spin for an infinite chain 
initialized in a symmetric Ramsey state with either no 
phase shift or a phase shift of A <p = ir between neighbor¬ 
ing atoms. The lattice constants have been chosen to be 
approximately Ao/2 as would be typical [5]. We refrain 
from choosing exactly Ao/2 to avoid the l/£ divergence. 
We observe that the dipoles’ lifetimes vary strongly, com¬ 
paring the subradiant behavior (red) where the collective 
dipole lives much longer than the natural lifetime of the 
atom to the superradiant (green) regime where the ex¬ 
citation vanishes very quickly. Additionally, to highlight 
the validity of the mean field approach, we add the results 
of a second order expansion simulation. Corresponding 
results for a full Ramsey sequence are shown in the ap¬ 
pendix D. 

Conclusions - In densely filled optical lattices dipole- 
dipole interaction and collective decay significantly 
change the evolution of an induced collective dipole. Due 
to the long range-nature of the coupling, sizable shifts ap¬ 
pear even for long lived clock states despite their minute 
dipole moment, which limits accuracy and precision of 
Ramsey spectroscopy. Shifts and dephasing in large sys¬ 
tems strongly depend on the dimensionality and geom¬ 
etry of the lattice, exhibiting resonant enhancements at 
particular lattice constants. While at current operating 
densities for Strontium [12-14] these shifts are smaller 
than other technical imperfections, they constitute in¬ 
herent fundamental perturbations even in perfectly filled 
lattice clocks. 

We have identified optimal operation geometries, 
which combine a negligible effective shift with a strong 
suppression of decay. In particular, for a ID lattice with 
a tailored excitation angle and for a 2D hexagonal lat¬ 
tice favorable operation parameters for future generation 
clock setups were found. In this sense it seems possible 
to implement a high density dark exciton based atomic 
clock geometry with shifts many orders of magnitude be- 
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FIG. 3. (Color online) Distance dependence of the effective quantities fi cff and r efi for an infinite equidistant chain, a square 
lattice and a hexagonal lattice (dashed black) compared to their not yet converged finite counterparts of 10, 4 * 10 4 and 10° 
particles respectively (solid red). Again, we find divergences at integral d/ Ao owing to the l/£-terms in F(£) and G(£)- In the 
2D configurations F eff plateaus at —1 for d < Ao, suggesting that this parameter range will be the most favorable for clock 
setups as decay is strongly suppressed. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 


a /Ao 

FIG. 4. (Color online) Effective interactions fl eff and r e£I for 
an infinite chain with spacing a where the spins are initially 
prepared with phase shift A <j> between neighboring spins. The 
dashed lines indicates parameters with f! cff = 0 optimal for 
an optical clock. 


low a single Hz and almost unlimited exciton life times. 
In 3D the interactions are particularly sensitive to the 
lattice constant and boundary effects, which dominate 



012345012345 
yt yt 

FIG. 5. (Color online) Three different examples for the time 
evolution of the spin expectation values for a chain with spac¬ 
ing d where initially all spins are prepared in a coherent su¬ 
perposition of ground and excited state with a phase shift of 
A <f>. The parameters used are d = 0.792Ao with A0 = 0 (red 
triangles), where D cff = 0 and r cfI is nearly optimal, as well 
as d = 0.49Ao (green squares) and d = 0.51Ao (blue circles), 
both with A0 = 7r which are close to a r efI discontinuity. The 
solid lines correspond to a solution of a second order cumulant 
expansion model with 200 particles and demonstrate a very 
good agreement with the infinite mean field description. 


even for billions of particles rendering such setups very 
challenging. 

While for most considerations we have focused on the 
case of clock transitions, the same physics is present 
in a more prominent and experimentally easier observ¬ 
able form for broader transitions. Optimizing geome¬ 
tries will also be relevant for devices such as superradiant 
lasers [32, 33] or lattice based optical memories. 
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Appendix A: Derivation of the mean field equations 

Starting from the full multiparticle density operator p our master equation allows to obtain the following equations 
for the individual spin expectation values immediately 


E x > = E 

^kj * 

b 

b 

\l (°%) + 

\ E 

(Ala) 

rd¥=k 



\iK) + 

rdAk 


A) = £ 

^kj * 


\ E T kj(<rj<r Z k) 

(Alb) 

rdAk 




rdAk 


y k ) = E 

1 


- <-K>) 

- 7(1 + K)) - ^ E Tk J (i a j a k) + ( a j a k)) ■ 

(Ale) 


jd^k jd^k 


Assuming a spatially separable state p = pk leads to the lowest order mean field equations used in the letter. 


Appendix B: Mean field equations with tailored excitation phase 

At zero temperature the ground state p = (£) k (\g){(g\) k is separable and in an idealized standard Ramsey procedure 
the first pulse would create a product state of equal weighted superpositions p = 1/2 ((|</) + |e))((e| + (g |)) fc . This 

is the generic initial state we use in our work to study dipole-dipole interaction. In fact, this state possesses the 
maximal collective dipole moment and therefore shows strong interactions. 

Of course, in any real setup this preparation step is not perfect as interactions are present during the excitation 
pulse and the excitation laser carries an intensity and phase gradient. Some of the errors can be corrected in improved 
excitation schemes[21, 22]. However, particularly in extended systems a phase gradient is hard to avoid and will 
strongly influence the system dynamics. Fortunately, one can show, that a known phase gradient will not complicate 
the calculations too much. If we allow for the individual atomic states to bare a spatially dependent phase of A<f> 
on the excited state, i.e. \ip k ) = ^ (Iff) + exp(*(/>*,)|e)), we can absorb this into our equations efficiently. Using the 
abbreviations Q c kj s = fl k j cos((j) k — (f>j) and fl s k ^ = fl k j sin(<j) k — </>j) we obtain the following modified equations of 
motion 


j t (aZ) = + E ^ 2 7( ^> + /> E -\Y, W> ( Bla ) 

jd¥=k j;j^k jd^k 

|k) = - e + E - \^i) + 2 E + \ E nr^^k) (Bib) 

JU^k jd^k jd^k jd¥=k 

j t wi) = - E + E <*”*%)) 

rd^k jd^k 

7(1 +{<)) \ e nri^^k) + <*k» - 2 E - <*k»- (bic) 

jd^k jd^k 

We see that the following definitions prove to be very helpful 

n k° S = E ^kjCOs((j) k - <t>j) n k n = E ^kj sin(^ fe - <j>j) (B2a) 


T c fe os = E cos ^k - h) r fe n = E F kj Mtk - 4>j) (B2b) 

rd^k rd¥^k 

Again, if we consider highly symmetric configurations where Q f = and T f = T| and the rotated states are initially 
identical we can define the effective rotated quantities 

n eS = n cos - ^r sin (B3) 

peff = pcos , 2 ^sin (B4) 
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which lead to a closed set of simplified effective equations as well, i.e. 

jfi x ) = rr'VXO - \l(Z x ) + \r eS (a x )(a z ) (B5a) 

j t {* v ) = -& s (a*)(a*) - ^(a v ) + (B5b) 

j t {° z ) = - 7(1 + W z )) ^r eff (V > 2 + <^} 2 ) (B5c) 

Note that such a phase gradient tends to mix the real and imaginary part of the interaction terms. 


Appendix C: Effective quantities for cubic lattices in 3D 

In a cubic 3D lattice the number of neighbors at a given distance r grows approximately as r 2 . Hence, one cam 
expect a slower convergence with distance. This problem is increased as the number of emitters to be considered 
grows with the third power of the system size. In contrast to ID and 2D, together these two scalings prevent a 
convergence of the effective interaction parameters in the range of tractable lattices sizes of up to N = (10 4 ) 3 = 10 12 
sites. Anyway, this is beyond experimentally realistic atom numbers so that we have to live with finite size effects. 

In order to demonstrate the very slow convergence of the infinite range mean field model, we present some typical 
intermediate result for a 3D cubic lattice. In Fig. 6 we depict the effective coupling strengths D cff and r cff for the 
innermost two-level system in a cubic lattice of about 8 • 10 9 particles, i.e. 2000 particles in each direction. We obtain 
strong and very rapid oscillations of the shifts as a function of the lattice constant. Notice, that 1/r-contributions as 
discussed in the letter will show up for planar and cubic diagonal distances of y/2 ■ r and \/?> ■ r as well. Increasing 
the atom number further still leads to changes of this pattern, so no final conclusions about physical properties and 
the behavior of a 3D cubic lattice can be obtained. However, perturbations of up to an order of magnitude larger 
than the linewidth as well as strong finite size shifts can be expected. In this case it is difficult to suggest an optimal 
lattice constant for a clock setup, except for avoiding certain resonances and choosing a region of about d ~ 3A/4. 




FIG. 6. Effective quantities and rf ff as experienced for the innermost spin inside a cube consisting of 2001 x 2001 x 2001 
spins in a cubic lattice configuration depending on the lattice spacing d. Even for very small changes of the lattice spacing the 
mean net-effect of all other spins will change dramatically. 


Appendix D: Ramsey spectroscopy 

The effective coupling and decay parameters H cff ,T efI characterize the interaction induced perturbation of the 
individual spin dynamics. Consequently, they will alter the Ramsey signal by introducing shifts of the fringes and 
modifications of the maximally obtainable contrast. As the actual connection between the magnitude of these effective 
couplings and their quantitative effect on the signal is nontrivial, we demonstrate the alterations of the Ramsey signal 



























































in the following examples. Using the previously derived equations of motion, it is straight forward to simulate the 
results of an ideal Ramsey sequence. By starting with a 7r/2-pulse all spins are rotated into the x-direction of the 
Bloch sphere. For a time yf the system evolves according to the equations 

<^> = -A aK) + U ef V>(^> - + *r ef V)<^> (Dla) 

{av) = \ a (al) - Q,° s (a*)(a z ) - W v ) + *r efl V)( CT *} (Dlb) 

W’) = - 7(1 + W z )) - eff (V ) 2 + (ay) 2 ) (Die) 

where A a = uiq — ujl is the detuning between the probe laser and the atomic transition frequency. After this free 
evolution a second 7r/2-pulse is applied and, finally, the expectation value of a z is measured. For a given system 
characterized by the effective quantities fi c£f and r efI the result of this measurement depends on the waiting time as 
well as on the detuning A a . In Fig. 7 the outcome of this numerical experiment is shown for three different realistic 
sets of effective quantities. The decisive quantity for the accuracy with regards to atomic clocks is the shift of the 
fringes due to the dipole-dipole interaction which can be obtained by measuring the shift of the maxima of the Ramsey 
fringes. The shifts for the chosen examples are shown in Fig. 8. On the other hand the slope of the fringes at their 
roots is the determining factor for the best achievable experimental precision. The numerical results are shown in 
Fig. 9. As seen in Fig. 10 the maximal shifts depend on ff cff only, while the maximal slope at zero points is governed 
by r eff . For realistic values for the effective quantities this means the accuracy can be limited to 7 and the achievable 
precision can vary by a factor of 5. 
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FIG. 7. Simulated idealized Ramsey spectroscopy for different sets of effective cooperative interaction strengths, from left to 
right: subradiant case without shift (S7 cff = 0,r efi = —0.75), independent atom limit (D cff = 0, r cfi = 0) and superradiant 
case with shift (fF® = l,r cff = 1). The colors indicate the free evolution time with cyan representing a very short time and 
magenta meaning times up to 2.57 _1 . 





FIG. 8. Shifts of the maxima after the free evolution time ■yt for different sets of effective interaction strengths as above, from 
left to right: (D eff = 0, T eff = -0.75), (D eff = 0,r eff = 0) and (D eff = l,F eff = 1). 
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FIG. 9. Slope of the signal at the zero crossing of the first fringe after a free evolution time of 7 1 for different sets of effective 
interaction strengths as above, from left to right: = 0, T cff = —0.75), (f2 eff = 0, r eff = 0) and = 1, r eff = 1). 
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FIG. 10. (a) Shift of Ramsey fringes depending on the effective coupling fl cff after t = 157 1 
. The different lines represent different choices of F cff , which hardly influence the result. The fringe shifts follow the effective 
mean field dipole coupling S7 cfI almost linearly and thus can be read off from the figures in the main manuscript, (b) 
Maximally achievable slope at roots depending on T ofi . The result is independent of the choice of Note that a negative 
r cff improves the measurement precision beyond the independent atom value. 


























